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Abstract 

We address the problem of constructing large undirected circulant 
networks with given degree and diameter. First we discuss the theo¬ 
retical upper bounds and their asymptotics, and then we describe and 
implement a computer-based method to find large circulant graphs 
with given parameters. For several combinations of degree and di¬ 
ameter, our algorithm produces the largest known circulant graphs. 

We summarize our findings in a table, up to degree 15 and diameter 
10, and we perform a statistical analysis of this table, which can be 
useful for evaluating the performance of our methods, as well as other 
constructions in the future. 
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1 Introduction 

The Degree/Diameter Problem (or DDP, for short) is one of the main theoret¬ 
ical issues in network design. DDP asks for constructing the largest possible 
graph, in terms of the number of vertices, subject to an upper bound on the 
maximum degree A and an upper bound on the diameter D. Let iVA,D be 
the number of vertices of the largest graph with maximum degree A and 
diameter D; it is well known that 
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-A^a,d ^ M/\ £) — 1 + A + A(A — 1) + ■ ■ ■ + A(A — 1)^ ^ 

^ f 1 + A if A > 2 

y 2D+ 1 if A = 2 

The number Ma,d is called the Moore bound, and a graph of order Ma,d 
is called a Moore graph [32]. Moore graphs are rather scarce: they exist only 
for a few combinations of A and D. For D = 1 and A > 1, they are the 
complete graphs on A + 1 vertices. For D > 2 and A = 2, they are the 
cycles on 2D +1 vertices. For D = 2, Moore graphs exist for A = 2, 3, 7, and 
possibly 57. 

The Degree/Diameter Problem can also be formulated for digraphs, but 
in this paper we will only be concerned with the undirected version. On 
the other hand, we are interested in some particular versions of the problem, 
namely when the graphs are restricted to a certain class, such as the class 
of bipartite graphs P [71 [HI HSl [E] , planar graphs |2Tl EH] HH [15] SHI Si] , 
vertex-transitive graphs [271 SO], Cayley graphs P E2l EZl SQII16], Cayley 
graphs of abelian groups m, or circulant graphs [El [Ml [29] . In particular, 
in this paper we are concerned with circulant graphs, i.e. Cayley graphs of 
hnite cyclic groups. 

For most of the aforementioned graph classes there exist Moore-like upper 
bounds, which are usually smaller than the Moore bound for general graphs. 
Most of the research, both in the general version, as well as in these restricted 
versions of DDF, falls into one of two main categories: 

1. Proving the non-existence of graphs with order close or equal to the 
upper bound, or 

2. Giving constructions of large graphs, whose order approach the upper 
bounds as much as possible. 

Research in the second category has been substantially boosted by the 
compilation of record graph tables, containing the largest known graphs for 
several combinations of degree and diameter. These tables provide bench¬ 
marks to test construction methods and computer search algorithms. Thus, 
after the hrst compilation started by Cornelias in 1995 [1], other record graph 
tables have been collected for general graphs, bipartite graphs, planar graphs, 
and Cayley graphs [5]. 

The study of circulant graphs began in 1970 with Elspas and Turner [T3] . 
Coincidentally, it was also Elspas who had formulated the Degree/Diameter 
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Problem back in 1964 na, but apparently he missed the connection between 
both topics. 

Even though the abelian property of the underlying group prevents abelian 
Cayley graphs in general (and circulant graphs in particular) to grow as large 
as their non-abelian counterparts, these graphs have been widely used as 
topologies for computer networks and parallel computers, due to their other 
nice properties. Paraphrasing m- . the extra structure provided by the 
groups may provide compensating advantages ..., such as good routing algo¬ 
rithms, easy constructibility, and the ability to map common problems onto 
the architecture”. 

However, up to now there were no comprehensive tables for abelian Cayley 
graphs, or their subfamily of circulant graphs. This is sufficient motivation 
to start a collection of benchmarks in the class of circulant graphs, which 
can be used for comparison purposes in the future. Our goals in this paper 
are twofold: 

1. Design and implement a computer search algorithm to find large cir¬ 
culant graphs with small degree and diameter, and 

2. Start the compilation of a table containing the largest known circulant 
graphs for some small values of degree and diameter, including the 
results obtained with the aforementioned algorithm. 

Additionally, we explore the sharpness of the existing upper bounds for 
some particular combinations of degree and diameter. 


2 Definitions and basic facts 

An undirected circulant graph C (n; S) is a Cayley graph on the cyclic group 
with a symmetric connection set S (i.e. S = S~^). Since is abelian, 
we can switch to additive notation and rephrase the symmetry condition of 
the connection set as S' = —S. In order to simplify the notation, we drop 
the curly brackets {} in the specification of the set S in C{n; S). 

Being a Cayley graph, C(n; S) is vertex-transitive. The degree of C(u; S) 
is A = |S|, and its order is obviously n. A circulant graph can also be dehned 
as a graph of n vertices whose adjacency matrix is circulant na. 

Regarding the degree, we distinguish two cases: 

1. Even degree: A = 2t. In that case, S = {±si,..., is*}, where 1 < 
Si< ... <St<^. 
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2. Odd degree: A = 2t + 1. In that case, S = {±si,..., ±St, |}, where 
1 < Si < ... < St < |. It follows that odd degree is only possible when 
n is even. 

0(n; S) is connected if, and only if, gcd(n, Si,..., s*) = 1. If gcd(?7,, r) = 1, 
then C{n] S) is isomorphic to C{n\ rS), where mnltiplication is taken modnio 
n. In that case we say that the connection sets S and rS are mnltiplicatively 
related. It should be noted, however, that two circulant graphs may be 
isomorphic without their connection sets being mnltiplicatively related [T^ 

E7]. 

Now let be the number of vertices of the largest circulant graph 

with degree A and diameter D. It was proved in [3] that, if A = 2f, then 

= ( 2 ) 

This upper bound was later rediscovered by Muga [36]. The quantity 
F(t,D) also turns out to be an upper bound for the order of the 

largest Cayley graph over an abelian group, with degree A and diameter D 
m- It is quite surprising that no better general upper bound (yet) exists for 
circulant graphs, considering that they are a special case of abelian Cayley 
graphs. 

The numbers F(t, D) of Eq. |^are known as Delannoy numbers (sequence 
A008288 of [3H|), and they arise in a variety of combinatorial and geometric 
problems |12|. For example, they correspond to the volume of the ball of 
radius D/2 in the metric in t dimensions [H] EU H5] . 

Unaware of the Delannoy connection, Stanton and Cowan had already 
studied these numbers back in 1970 [H], as a generalization of the binomial 
coefficients, which satisfy the recurrence: 


F{t, D) = F{t — 1, D) + F(f, D — 1) + F{t — 1, D — 1), with 
F(f, 1) = 2t + 1, for f > 0. 


They gave several interesting formulas for these numbers, such as: 



(4) 


In the case of odd A (i.e. A = 2f + 1) we have the generator |, which 
is its own inverse. Figure provides a graphical example of that case, for 
A = 5 and D = 3. We have denoted the generators as a, b, c, where c = —c. 
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Figure 1: Tree representation of a maximal abelian Cayley graph of 
degree 5 and diameter 3. 


We can see in this example that the edges labeled with c duplicate every 
vertex, except those in the lowest level (level D). Therefore, an upper bound 
for the maximum number of vertices in this case is: 

F\t,D) = F{t,D) + F{t,D-l) 

From Equations and we can deduce that the numbers F'(t, D) 
satisfy the same recurrence: 

F'(t, D) = F'{t - 1, T>) + F'(t, D-l) + F\t - 1, D - 1), with 
F'(t,l) = 2t + 2, fort>0. 

Additionally, F{t,D) < F'{t,D) < F{t + 1,-D), i.e. the sequence of 
upper bounds obtained by fixing the diameter is monotonically increasing. 
Nevertheless, F(t, D) and F'{t, D) grow at slightly different rates, as we will 
see next. 

In order to hnd an explicit formula for F'{t,D), we proceed as in [3T] : 
Dehne the generating function Ad^z) = Multiplying both 

sides of Eq. by z* and summing over f > 1 we get 


(5) 


( 6 ) 


Ad{z) — A£)(0) — zAd{z) + {Ad^i^z) — A£)-i(0)) + zAb-i^z) (7) 
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whence 


^d{z) = YZr^Ao-iiz) 

With the aid of the boundary condition in Eq. [^we get 


( 8 ) 


Ad{,z) 


2(1 + 
(1 — 


(9) 


Now, A]:){z) is the product of (1 + 2:)'^ ^ t l/{l — z)^^^ = 

~ Si Then, the series of Ad{z) can be obtained as the 

convolution of the respective factor series, which gives us: 


F\t,D) 



( 10 ) 


Note that Eq. 


10 is the same as Eq. 2 of ra, A and D 

swapped]^ This symmetry arises as a consequence of the symmetry of Eq. 

In the same manner we can derive other formulas for F(f, D) from Eq. 2 of 



With the aid of the methods developed in |18] it can be shown that the 
numbers F{t, D) and F'{t, D) do not have a closed form, meaning that they 
cannot be represented as a linear combination of a fixed number of hyperge¬ 
ometric terms. However, we can obtain asymptotic formulas for them. From 

^More precisely, t takes the value p (which corresponds to [i2/2j in [SJ), and D takes 
the value k (which corresponds to A/2). 
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Corollary 2.1 of 
both F(t,D) and F'(t,D) grow asymptotically as 


31], and nsing the symmetry jnst discnssed, we get that 

(2t)^ 


D\ 


+ o{t 


D\ 


np 


The hrst few valnes of F{t,D) and F'{t,D) are collected in Table ^ 
to degree 15 and diameter 10. We can £t an exponential fnnction to 

these valnes, where 


/(A, D) = - 0.8015 + 0.5612A + 0.5433T> - 0.05421A2 + 0.1499AT) 

- 0.09784T)2 + 0.001574A^ - 0.001924A2T) (12) 

- 0.003822AT)2 + 0.004591T>^ 

with = 0.9986. 

Fignre|^ shows a contonr plot of /(A, 77). 


<S> Contour plot for logarithms of upper bounds 
• Logarithm of upper bound vs. Degree and Diameter 
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Figure 2: Contour plot of the function f{A,D), that approximates 
the logarithms of the upper bounds F{t, D) and F'{t, D). 

If we compute the differences between the values predicted by f{A,D) 
and the logarithms of the actual upper bounds F(t,D) and F'(t,D), we 
see that the numbers F{t,D) are slightly above f{A,D), while F'{t,D) are 
slightly below /(A, D), as shown in Figure]^ Therefore, it might be useful to 
compute different approximations for the values with even and odd degree, 
respectively: 
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/(A, D)=- 0.5686 + 0.4266A + 0.6491D - 0.03746A2 

+ 0.1366AL) - 0.1009L)2 ^ o.00098A3 (13) 

- 0.001492A2L) - 0.003578AZ12 ^ 0.004598D^ 

(with = 0.9999), and 

/(A, D) = - 0.6153 + 0.5063A + 0.4993D - 0.04712A2 

+ 0.1507AZ1 - 0.09576L)2 + 0.001282A3 (14) 

- 0.001776A2Z1 - 0.003974AZ12 ^ 0.004584/1^ 


with = 0.997. 



Figure 3: Predicted (logarithmic) upper bound minus actual (loga¬ 
rithmic) upper bound. Light entries correspond to nega¬ 
tive differences. Darker colors correspond to positive dif¬ 
ferences. 

If we normalize the difference at each (A, D) by the logarithm of the 
actnal npper bonnd at {A,D), the disparity between even and odd rows is 
still visible, bnt it fades ont as A increases, as shown in Fignrej^ 

Circnlant graphs of the form C{n-, ±1, ±S 2 ,..., ±St) are called multi-loop 
graphs. In particnlar, for t = 2 and t = 3 they are called double-loop graphs 






Figure 4: Normalized difference between the predicted (logarithmic) 
upper bound minus actual (logarithmic) upper bound. 
Light entries correspond to negative differences. Darker 
colors correspond to positive differences. 
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and triple-loop graphs, respectively. According to |35], the maximnm order 
of a triple-loop network C(n; ± 1 , ±S 2 , ± 53 ) is: 


T< 


6,D 


§A)3 + mjj2 + 2D + 1 

= <; 32[Sj=+48[fJ^ + 30Lej+7 
32L|j=‘ + 80[|j'^ + 70LfJ+21 


if D = 0 mod 3 
ii D = 1 mod 3 
if D = 2 mod 3 


Many of the graphs fonnd by ns here are mnlti-loop. 


(15) 


3 A basic search procedure 

Let n, A and D be given. We now describe a systematic search procednre 
that attempts to hnd a circnlant graph with the given parameters. Onr search 
space consists of all connection sets S. The basic idea is to consider a tree-like 
organization of the search space, and then perform a Depth-First Search with 
backtracking on the tree, while prnning signihcant portions of it. To begin 
with, onr search space will be restricted to connection sets with si = 1. This 
will rednce the search space withont losing too many potential solntions, and 
as a bonus, it will relieve us from the burden of testing connectedness. 

For even (resp. odd) A, a node at depth t in the search tree corresponds 
to a circnlant graph of order n and degree 2t (resp. 2f -|- 1). For the sake of 
brevity, we shall make no distinction between a node and its corresponding 
graph. The root of the search tree is the circnlant graph C(? 7 ,;±l}) (resp. 
C(? 7 ,; {± 1 , n/ 2 )). Let C{n;S) be a node of degree smaller than A, and let 
m be the largest element in S such that 1 < m < n/2. Then C{n] S) has 
offspring C{n] S U {±( 7 }) for every m + k < g < {n — k)/2 (resp. m + k < 
g < n/2 — k), where k is a. constant standing for the maximum difference 
allowed between two elements in S. Nodes of degree A, appearing at depth 
(A — l)/2 (resp. (A — 3)/2), are the leaves of the tree. 

In the general Degree/Diameter Problem, it is well-known that large 
graphs tend to have large girth with respect to their diameter. This ob¬ 
servation does not apply to circnlant graphs, as their girth is always at most 

4 when A > 3. Thus, we have generalized the girth criterion as follows: 

For a circnlant graph C{n; S) and a length d, dehne {0 < w < N) as 
the set of paths of length d from a hxed vertex v (say 0) to the vertex w. Paths 
in are meant to be disjoint, with no repeated vertices (except when w = 
v), and they are different up to commutation (that is, P = Qx{xa){xab)R 
and P' = Qx{xb){xha)R = Qx{xb){xab)R are considered to be the same 
path). 
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Whether the offspring of a given node C{n] S) are explored or not depends 
on a function on the cardinals of the sets Basically, for each i and d we 
hx constants cf (for 0 < i < d < D) which stand for the maximum number 
of paths allowed in whenever w is at distance i < d from the vertex 0. 
That is, for every vertex w in the graph, if w is at distance i < D from the 
vertex 0, we require that \V!^\ < cf for every d with i < d < D. 

Tuning the constants cf provides great flexibility on the number of nodes 
explored, and proved to be very effective in practice. By defining small 
enough values for the constants cf it is possible to prune a significant portion 
of the search tree, especially as A increases. On the other hand, the constants 
cf should not be too small, so as not to miss too many potential solutions. 

Therefore, a careful choice of the constants cf is of paramount importance. 
In order to tune them properly, some degree of experimentation and intuition 
was involved. Initially, we tested many of the largest known circulant graphs 
for the values of \Vfj \. Most of these graphs were already known to be optimal. 

The parameter k, introduced to control the size of the search space, and 
thus the number of nodes visited by the algorithm, presented similar trade¬ 
offs. Although certainly useful, experimenting with different values for the 
parameter k provided rather poor improvements, compared to the cut-offs 
achievable via pruning (unless the ability to hnd existing solutions was seri¬ 
ously compromised by choosing k to be too large). 

Algorithm [T] formalizes the above ideas. Parameters n, A and D are 
global. The input to the algorithm is the connection set S. Initially, the 
algorithm is called with S = {±1} (resp. S = {±l,n/2}). Every time that 
a new generator is added, the algorithm calls itself recursively with the new 
connection set S'. If a circulant graph C{n-, S) with de desired parameters is 
found, the algorithm will print it out. Otherwise, the algorithm terminates 
without producing any output. 

Although the algorithm is presented in a recursive fashion, it was actu¬ 
ally implemented using a stack for optimisation reasons. We omit actual 
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implementation details here. 
Algorithm 1: SEARCH 


Input : A set of generators S. 

Output: A circulant graph C(n; S) with degree A and diameter D. 


3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 


1 m:= Largest element in S'; 

2 for q:= m + k to in — k)/2 do 

success := true; 

for each g-path P in C(n; S'), with P e do 
i:= Distance from vertex 0 to tn in C{n] S')] 


if d > i and \V!^\ > cf then 
success := false; 

break; 
end 

if d = i and \Vi,\ > cl for some j (i < j < D) then 
success := false; 

break; 
end 
end 

if success then 
if |S"| < A then 

I Call SEARCH(^'); 

end 
else 

if diameter of C{n] S') is D then 
Print C(n; S'); 

end 
end 
end 


26 end 


With the aid of Algorithm we have been able to hnd 15 new large 
circulant graphs of small diameter (ranging from 3 to 5), and degrees between 
8 and 15. Additionally, we found two circulant graphs, which although not 
being record graphs, will be of interest later, since they will be used in the 
construction of larger graphs. The connection sets for all these graphs are 
listed in Table Four of them were found independently by Lewis in [25] , 
and proved to be optimal. These four graphs, as well as other graphs in 
Table up to degree 13, were also reported in [TH|. By exhaustive search we 
also prove the optimality of the graph with degree 11, diameter 3, and order 
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210, as well as the optimality of the previously known graphs of diameter 3 
and degrees 9 and 10. 


A 

D 

Order / Connection Set 

8 

3 

104 {±1, ±16, ±20, ±27} 



4 

248 {±1, ±61, ±72, ±76} 



5 

528 {±1, ±89, ±156, ±162} 




511 {±1, ±5, ±70, ±96} 



6 

967 {±1, ±7, ±132, ±182} 


9 

4 

320 {±1, ±15, ±25, ±83,160} 


10 

4 

457 {±1, ±20, ±130, ±147, ±191} 

★ 


5 

1099 {±1, ±53, ±207, ±272, ±536} 

L 

11 

3 

210 {±1, ±49, ±59, ±84, ±89,105} 

★ 


4 

576 {±1, ±9, ±75, ±155, ±179, 288} 

★ 


5 

1380 {±1, ±33, ±173, ±387, ±663,690} 

★ 

12 

3 

275 {±1, ±16, ±19, ±29, ±86, ±110} 

★ 


4 

761 {±1, ±12, ±184, ±235, ±334, ±362} 

★ 


5 

1800 {±1, ±30, ±64, ±384, ±761, ±841} 

★ 

13 

3 

312 {±1, ±14, ±74, ±77, ±130, ±138,156} 

★ 


4 

920 {±1, ±11, ±38, ±176, ±232, ±376,460} 

★ 

14 

3 

381 {±1, ±11, ±103, ±120, ±155, ±161, ±187} 

★ 

15 

3 

448 {±1, ±10, ±127, ±150, ±176, ±189, ±217, 224} 

★ 


Table 1: Connection sets for graphs found by Algorithm labelled 
with a Graphs labelled with an L were found indepen¬ 
dently by Lewis [25], and proved to be optimal. Graphs 
labelled with a _ symbol are sub-optimal. The graph with 
order 1099 was also found by Lewis [26j . 


Unfortunately, for larger values of A and D, the execution time of Al¬ 
gorithm becomes prohibitive with the computing power at our disposal. 
Thus, in order to complete our table of record circulant graphs up to A = 15 
and H = 10, we must use a combination of Algorithm [T] with Cartesian 
product, as described in the next section. 

4 Combining basic search with Cartesian prod¬ 
uct 

Prior to our computer-based search, the main sources of large circulant 
graphs were: 
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1. For t = 2 an optimal circulant graph (^(n; ±si, ± 52 ), is achieved for 

Si = - 1 - 1)J and S 2 = Si + 1 [Tl 1^ |33]. 

2. Monakhov and Monakhova nsed an evolntionary algorithm to hnd 
dense families of nndirected circnlant graphs. In particnlar, with the 
aid of this algorithm they fonnd some families of large triple-loop graphs 

H- 

3. For larger degrees we have the constrnction (^(n; ±1, ±s,..., ±s^“^), 

where n = s^, and s is an odd integer, which yields good circnlant 
graphs of diameter |(s — 1) = | |17j. 

4. Applying the methods described in |TT] for abelian Cayley graphs, 
Charles Delorme has recently obtained several circnlant graphs of large 
order uni. 

5. With the aid of compnter search, Lewis fonnd several optimal graphs 
with degrees 8 and 9 [25]. 

6. With the aid of Algorithm]^ we fonnd some new record graphs, np to 
degree 13 [T8] . 

Onr new compnter search method is based on the combination of Algo¬ 
rithm above with Cartesian prodnct of graphs. If n and m are relatively 
prime, then x Given the Cayley graphs C{n] Si) and G(m; S' 2 ), 

the Cayley graph C{nm] S) can be obtained as C{n] S'!) □ C{m] S' 2 ), where 
S = mS'i U nS' 2 , and □ represents the Cartesian prodnct of graphs. In this 
constrnction, the degree of C[nm] S) is the snm of the degrees of the fac¬ 
tor graphs C{n]Si) and C{m;S 2 ), and the same holds for the diameter of 
C{nm;S). Figure depicts the Cartesian product C(4;l) □ (7(3;!), and 
Figure shows this same graph with a circulant layout. Note that the gen¬ 
erator 3 of Z 12 corresponds to (1, 0) in the -external- direct product of Z 4 
and Z 3 , while the generator 4 corresponds to (0,1). 

Now we just have to look for appropriate combinations of the largest 
known graphs in order to £11 higher entries of the table. Normally we would 
be interested in using the largest possible factor graphs, but in order to apply 
the aforementioned construction we may need to resort to sub-optimal factor 
graphs, since the optimal factors may have orders which are not relatively 
prime, hence yielding a non-circulant graph. 

Algorithm formalizes the method. The auxiliary function prime{G) 
determines whether G is the Cartesian product of two circulant graphs whith 
relative prime orders. With the aid of Algorithm we found 41 new record 
circulant graphs, listed in Table 
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Figure 5: Cartesian product (7(4; ±1) □ (7(3; ±1) 

0 
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Figure 6: The graph (7(12; ±3, ±4) 
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Algorithm 2: COMBINED SEARCH AND CARTESIAN PRODUCT 


Input : A degree A > 4 and a diameter D > 10. 

Output: Largest graph C7(nm; S) = C{n; Si)OC{m] S 2 ) with degree A 
and diameter D; where C{n] Si) and (7(m; S 2 ) are known 
circulant graphs, or are found by Algorithm [T] 

1 max := 0; 

2 A:= Empty array of graphs; 

3 for i:= 2 to [A/2J do 

4 

5 

6 

7 

8 
9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 


for j:= 1 to [D/2J do 

Gi := Largest known circulant of deg. i and diam. j; 

G 2 := Largest known graph of deg. A — i and diam. D — j; 
if IG 1 IIG 2 I > max then 

Insert |Gi|Q|G 2 | in A so as to keep A sorted; 
if gcd{\Gi\, IG 2 I) = 1 then 
max := |Gi||G 2 |; 

Remove from A all the elements after |Gi|D|G 2 |; 

end 
end 
end 
end 

while lenqth(A) > 0 and not prime(A\0]) do 
G := |Gi|n|G2| :=^[0]; 

Remove G from A-, 
n' := |Gi| - 1; 
while n'\G 2 \ > max do 
if gcd{n\ m) = 1 then 

G[ := Circulant of order n', same deg. and diam. as Gi, 
if G[ was found then break; 
n' := n’ — 1; 

end 
end 

m’ ■= IG 2 I - 1; 
while |Gi|m' > max do 
if gcd{n, m') = 1 then 

G '2 := Circulant of order m', same deg. and diam. as G 2 ] 
if G 2 was found then break; 
m' := m' — 1; 

end 
end 

G' := n'm > nm' ? G'^UGi : GiUG' 2 ] 

Insert G' in A so as to keep A sorted; 
max := |G'|; 

Remove from A all the elements after G'; 


end 

40 if length{A) > 0 then Print ^[0]; 
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Gi 

G2 

A 

D 

Order 

A 

D 

Order 

A 

D 

Order 

10 

6 

1533 

2 

1 

3 

8 

5 

511 


7 

2 925 

4 

3 

25 

6 

4 

117 


8 

5136 

2 

1 

3 

8 

7 

1712 


9 

8 560 

2 

2 

5 

8 

7 

1712 


10 

13 840 

2 

2 

5 

8 

8 

2 768 

11 

6 

2100 

2 

1 

3 

9 

5 

700 


7 

4 088 

3 

2 

8 

8 

5 

511 


8 

7 736 

3 

2 

8 

8 

6 

984 


9 

13400 

4 

3 

25 

7 

6 

536 


10 

21976 

4 

4 

41 

7 

6 

536 

12 

6 

3 297 

2 

1 

3 

10 

5 

1099 


7 

6 864 

4 

2 

13 

8 

5 

528 


8 

13 200 

4 

3 

25 

8 

5 

528 


9 

24 600 

4 

3 

25 

8 

6 

984 


10 

42 800 

4 

3 

25 

8 

7 

1712 

13 

5 

1828 

3 

1 

4 

10 

4 

457 


6 

4 396 

3 

1 

4 

10 

5 

1099 


7 

9100 

4 

2 

13 

9 

5 

700 


8 

18 720 

6 

4 

117 

7 

4 

160 


9 

36 036 

6 

4 

117 

7 

5 

308 


10 

63 700 

4 

3 

25 

9 

7 

2 548 

14 

4 

825 

2 

1 

3 

12 

3 

275 


5 

2 285 

4 

1 

5 

10 

4 

457 


6 

5 941 

4 

2 

13 

10 

4 

457 


7 

14 287 

4 

2 

13 

10 

5 

1099 


8 

29 016 

6 

4 

117 

8 

4 

248 


9 

59 787 

6 

4 

117 

8 

5 

511 


10 

113139 

6 

4 

117 

8 

6 

967 

15 

4 

1100 

3 

1 

4 

12 

3 

275 


5 

3 044 

3 

1 

4 

12 

4 

761 


6 

7 524 

7 

3 

76 

8 

3 

99 


7 

17940 

4 

2 

13 

11 

5 

1380 


8 

39 564 

5 

3 

36 

10 

5 

1099 


9 

81900 

6 

4 

117 

9 

5 

700 


10 

154 720 

7 

4 

160 

8 

6 

967 

16 

5 

3 805 

4 

1 

5 

12 

4 

761 


6 

10 296 

8 

3 

99 

8 

3 

104 


7 

25135 

6 

3 

55 

10 

4 

457 


8 

60445 

6 

3 

55 

10 

5 

1099 


9 

128 583 

6 

4 

117 

10 

5 

1099 


10 

269 808 

8 

ft 

511 

8 

5 

528 


Table 2: Large circulant graphs obtained by the Cartesian product 
of two circulant graphs. 




5 The largest known circulant graphs 

Table 1^ summarizes our current knowledge about the largest known circulant 
graphs, including the ones obtained with the aid of our methods. Each entry 
of the table contains four elements: the order of the largest circulant graph 
known to-date (upper left), a reference to its source (upper right), the best- 
known upper bound (lower left), and the percentage of the upper bound 
achieved (lower right). This table can also be found in [B]. 

The actual graphs obtained by Algorithms and are given in tables 
and The tables also indicate which of these graphs are optimal multi-loop 
graphs. Optimality has been verihed for some small values of A and D, by 
performing a systematic search through all multi-loop graphs with the given 
parameters A and D. 

Table can be well approximated by a bi-exponential function, such as 
where 

/(A, D)=- 1.912 + 0.8845A + 0.8014/1 - 0.08641A2 

+ 0.1152A/1 - 0.1214/12 ^ 0.002799A3 (16) 

+ 0.002799A2/1 - 0.001046A/12 + 0.005036/l^ 

with i/2 = 0.997. 

Figure shows the approximation of the bi-cubic polynomial /(A,/I) to 
the logarithms of the graph orders. 

Fitting a function to the table provides insight in several ways. Obvi¬ 
ously, it gives an idea on how the table grows, and how effective the different 
construction methods are. Second, from this approximation we can hgure 
out which are the table entries where an improvement is more likely to oc¬ 
cur in the near future, because their current value lies beneath the value 
predicted by the function Figure summarizes this analysis in a 

graphical way. Entry (A, D) contains a color representation of the predicted 
value f{A,D) minus the actual value of the logarithmic (A,/I) table. The 
lighter the color, the better the actual value is, in comparison with the value 
predicted by f{A,D). Dark entries correspond to ‘bad’ values, which lie 
below their prediction. 

There are several interesting things that pop out immediatlely: First of 
all, odd rows tend to be darker than even ones, especially for lower degrees. 
However, this should be no surprise to us, since the upper bounds show the 
same pattern, as we saw in Section When we look at the percentages 
(next), we see that even rows actually do not fare better than odd ones. 

The best opportunities for improvement seem to lie in rows 10 and 11, 
from diameters 6 to 10, and also in the small black region of rows 14,15, 
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Table 3: Orders of the largest known circulant graphs. 



















A 

D 

Circulant graph 

7 

6 

536{±1,±231,±239,268} * 

7 

7 

828{±1,±9,±91,414} * 

7 

8 

1232{±1,±11,±111,616} * 

7 

9 

1764{±1, ±803, ±815,882} * 

7 

10 

2392{±1,±13,±183,1196} * 

8 

3 

104{±1,±16,±20,±27} ★ 

8 

4 

248{±1,±61,±72,±76} ★ 

8 

5 

528{±1,±89,±156,±162} 

8 

6 

967{±1,±7,±132,±182} 

8 

7 

1545{±1, ±170, ±178, ±468} 

9 

4 

320{±1, ±15, ±25, ±83,160} * 

9 

5 

684{±1, ±111, ±145, ±279, 342} 

9 

6 

1284{±1, ±36, ±163, ±342, 642} 

9 

7 

2340{±1, ±149, ±157, ±645,1170} 

10 

4 

457{±1, ±20, ±130, ±147, ±191} 

10 

5 

1099{±1, ±53, ±207, ±272, ±536} 

10 

6 

1533{±3, ±15, ±210, ±288, ±511} 

10 

7 

2925{±25, ±351, ±400, ±468, ±550} 

10 

8 

5136{±3, ±645, ±1712, ±1824, ±1848} 

10 

9 

8560{±5, ±1075, ±1712, ±3040, ±3080} 

10 

10 

13840{±5, ±1032, ±2768, ±5360, ±5400} 

11 

3 

210{±1, ±49, ±59, ±84, ±89,105} ★ 

11 

4 

576{±1, ±9, ±75, ±155, ±179, 288} 

11 

5 

1380{±1, ±33, ±173, ±387, ±663,690} 

11 

6 

2100{±3, ±15, ±591, ±669, ±700,1050} 

11 

7 

4088{±8, ±40, ±511, ±560, ±768, 2044} 

11 

8 

7736{±8, ±56, ±967, ±1056, ±1456,3868} 

11 

9 

13400{±25, ±1608, ±2144, ±5775, ±5975, 6700} 

11 

10 

21976{±41, ±2144, ±2680, ±9471, ±9799,10988} 

12 

3 

275{±1, ±16, ±19, ±29, ±86, ±110} 

12 

4 

761{±1, ±12, ±184, ±235, ±334, ±362} 

12 

5 

1800{±1, ±30, ±64, ±384, ±761, ±841} 

12 

6 

3297{±3, ±159, ±621, ±816, ±1099, ±1608} 

12 

7 

6864{±13, ±1056, ±1157, ±1584, ±2028, ±2106} 

12 

8 

13200{±25, ±1584, ±2112, ±2225, ±3900, ±4050} 

12 

9 

24600{±25, ±2952, ±3936, ±4075, ±8700, ±8850} 

12 

10 

42800{±25, ±5136, ±5375, ±6848, ±15200, ±15400} 


Table 4: New circulant networks obtained by Algorithms and 


The graphs marked with * are optimal. 
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A 

D 

Circulant graph 

13 

2 

80{±1, ±3, ±9, ±20, ±25, ±33,40} 

13 

3 

312{±1, ±14, ±74, ±77, ±130, ±138,156} 

13 

4 

920{±1, ±11, ±38, ±176, ±232, ±376,460} 

13 

5 

1828{±4, ±80, ±457, ±520, ±588, ±764, 914} 

13 

6 

4396{±4, ±212, ±828, ±1088, ±1099, ±2144, 2198} 

13 

7 

9100{±13, ±65, ±1400, ±2100, ±2561, ±2899,4550} 

13 

8 

18720{±117, ±160, ±585, ±2560, ±3520, ±3627,9360} 

13 

9 

36036{±117, ±308, ±819, ±4928, ±5031, ±6776,18018} 

13 

10 

63700{±25, ±175, ±7644, ±10192, ±13025, ±14275,31850} 

14 

2 

90{±1, ±4, ±10, ±17, ±26, ±29, ±41} ★ 

14 

3 

381{±1, ±11, ±103, ±120, ±155, ±161, ±187} * 

14 

4 

825{±3, ±48, ±57, ±87, ±258, ±275, ±330} * 

14 

5 

2285{±5, ±100, ±457, ±650, ±735, ±914, ±955} * 

14 

6 

5941{±13, ±260, ±914, ±1371, ±1690, ±1911, ±2483} * 

14 

7 

14287{±13, ±689, ±2198, ±2691, ±3297, ±3536, ±6968, } * 

14 

8 

29016{±117, ±248, ±3968, ±5456, ±7137, ±8424, ±8892} * 

14 

9 

54120{±55, ±984, ±4920, ±8965, ±19140, ±19470, ±20664} 

14 

10 

113139{±117, ±819, ±967, ±15444, ±15472, ±21274, ±21294} 

15 

2 

96{±1, ±3, ±5, ±11, ±24, ±31, ±39,48} 

15 

3 

448{±1, ±10, ±127, ±150, ±176, ±189, ±217, 224} * 

15 

4 

1100{±4, ±64, ±76, ±116, ±275, ±344, ±440, 550} 

15 

5 

2880{±5, ±45, ±375, ±576, ±775, ±895, ±1152,1440} 

15 

6 

7488{±13, ±117, ±975, ±1152, ±1728, ±2015, ±2327, 3744} 

15 

7 

17584{±16, ±848, ±1099, ±3297, ±3312, ±4352, ±8576, 8792} 

15 

8 

39564{±36, ±1099, ±1908, ±5495, ±7452, ±9792, ±19296,19782} 

15 

9 

81900{±117, ±585, ±700, ±11200, ±15400, ±23049, ±26091,40950} 

15 

10 

154720{±160, ±967, ±1120, ±4835, ±21120, ±29120, ±29977, 77360} 


Table 5: New circulant networks obtained by Algorithms and 


The graphs marked with * are optimal. 
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Degree 


Contour plot of cubic fit 
• Log of order vs. Degree and Diameter 



5 6 7 

Diameter 


Figure 7: Contour plot of the bi-cubic polynomial f{A,D), com¬ 
pared to the logarithmic Degree/Diameter table 


Diameter 



Figure 8: Predicted (logarithmic) table values minus actual (loga¬ 
rithmic) values. Light entries correspond to negative dif¬ 
ferences. Darker colors correspond to positive differences. 
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columns 4, 5. On the other hand we get a light triangular area in the low 
diameters, in rows 10 — 15. That is where we have been able to perform a 
more detailed search with Algorithm 

Another important analysis tool is given by the table of percentages, 
which indicates for each A and D the percentage of the upper bound that 
has been achieved by the largest known graph. Figure shows the table of 
percentages as a contour plot. 



Figure 9: Table of percentages shown as a contour plot 

It is not unreasonable to assume that the percentages should behave 
smoothly, therefore we can subject the table of percentages to the same 
analysis that we have used for the table of upper bounds and the table of 
graph orders. In this case, a bi-cubic polynomial is not flexible enough to 
approximate the table of percentages, therefore, we have htted the bi-quartic 
polynomial 

g(A, D) = 98.77 - 13.16A + 17.23T> + 1.214A2 + 4.512AT) 

- Q.777D‘^ + 0.03356A^ - 0.89A^T) + 9.2892AD‘^ 

+ 0.6634T)^ - 0.003327A^ + 0.02822A^T) ^ ’ 

+ 0.00937A2T)2 _ o.01317AT)3 - 0.02335T>^ 

to the table of percentages. The contour plot of g(A, D) is given in Figure 

uni 
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Quartic fit to percentages 
• Percentage vs. Degree and Diameter 



Diameter 

Figure 10: Contour plot of q{A,D), which approximates the table 
of percentages 


Now, for each A and D we can constrnct the difference between the 


predicted valne q{A,D) and the actnal percentage. Fignre 11 shows these 
differences, normalized by the actnal percentages. Again, the lighter the 
color, the better the graph is, in comparison with its expected percentage. 


Fignre 11 conhrms the main conclnsions drawn from the analysis of Fignre 


6 Open problems 

The apparent simplicity of circnlant graphs is deceiving, and they are a 
reservoir of interesting open problems. The main theoretical qnestion here 
has to do with hnding sharper npper bonnds for circnlant graphs. The bonnds 
F(t, D) and F'{t, D) given by Eqnations ii [l0| and are clearly not 
sharp. All graphs in Table from degrees 5 to 8 are known to be optimal, 
and yet they lie below the npper bonnds. Some sharper lower bonnds have 


been fonnd for specihc degrees or diameters, snch as Eq. bnt no general 
sharper lower bonnd is known np to date. 

The npper bonnds F(f, D) and F'{t, D) have been obtained only by nsing 
the commntativity of the nnderlying gronp. However, circnlant graphs are 
not even optimal in the class of abelian Cayley graphs (see m and also HI), 
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Diameter 

23456789 10 




0 

- 0.05 


- 0.1 


- 0.15 


Figure 11; Normalized differences between the predicted percent¬ 
ages and the actual percentages. Light entries corre¬ 
spond to negative differences. Darker colors correspond 
to positive differences. 
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despite some papers claiming so m- 

As for the lower bounds, there are few constructions for A > 8. Our 
Cartesian product construction may clearly be generalized to obtain a general 
lower bound, which may then be compared analytically with other existing 
constructions. 

Regarding computer-based search, this problem poses several challenges. 
The search space grows very quickly with A and D, making a systematic 
search unfeasible. This problem could be alleviated by combining the search 
procedure with some isomorph rejection sieve, or some powerful heuristics. 
Additionally, other computer-based methods in the style of [M] could be 
devised, to obtain families of large circulant graphs with low diameter. 

In this paper we have only addressed the construction of undirected cir- 
culants, but in principle, the methods described here can be extended to 
directed, or even mixed graphs. 
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